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We investigate numerically the dynamically generated plastic deformations of a 3D 
vortex lattice (VL) driven through a disorder potential with isolated, strong pinning 
centers (point-like or extended along the field direction). We find that the VL 
exhibits a very peculiar dynamical behavior in the plastic flow regime, in particular, 
topological excitations consisting of three or four entangled vortices are formed. We 
determine the critical current density j c and the activation energy for depinning U c 
in the presence of a finite density of strong pinning centers. 

PACS numbers: 74.60.-w, 74.60.Ge 
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Recently, the dynamical behavior of a VL in a disordered type-II superconductor has 
attracted much interest . The intriguing dynamics of the VL originates from the compe- 
tition between the vortex - defect and the vortex - vortex interaction, leading to a threshold 
behavior If the pinning potential dominates then the vortices are essentially stationary 
with a slow residual motion due to tunneling or thermal activation. In the opposite limit, 
when the pinning forces are weak compared with the driving force, the VL is only elastically 
deformed and flows coherently. In the important intermediate regime, when the pinning and 
driving force are comparable, the VL deforms plastically and an incoherent motion results 
HH . The breakdown of the weak collective pinning scenario |7]] with increasing strength of 
the pinning potential has been observed in numerical simulations ||. Furthermore, it has 
numerically been shown |||9| that the plastic flow of a 2D VL driven through a random 
pinning potential consists of channels of flowing vortices, i.e., rivers, and regions of pinned 
vortex lines, i.e., islands. This so-called channel flow behavior has been observed experimen- 



tally by means of Lorentz microscopy pTj|JTT|] . The various anomalies P,|12j (e.g., thermal 



instability at large currents, non-monotonicity of the I-V curves, peak effect) occurring in 
the plastic flow regime suggest that the nature of the dynamics is fundamentally dissimilar 
from the one in the elastic regime. In particular, the dynamically generated disorder is 
of general importance as it seems to be a basic constituent of the dynamics in disordered 
systems. 

In this paper we present for the first time a study on the dynamically produced defects 
in a 3D VL. We concentrate on a VL driven through a material with isolated strong pinning 
centers. It turns out that such strong pinning centers can entail the generation of localized 
topological excitations involving three or four entangled vortices. Moreover, the presence 
of strong pins leads to a preferred orientation of motion of the VL and affects the critical 
current density as well as the barrier for depinning. In the following we first introduce our 
model and then consider point-like as well as extended pinning centers. 

If the density of defects is low, the VL deformations produced by individual pinning 
centers can be studied without accounting for collective effects. The pinning centers con- 



sidered here are assumed to have a rod-like shape with a small lateral (coherence length 
~ £) but a variable longitudinal (from £ to a few lattice constants a® along the field direc- 
tion) extension. This kind of defects corresponds, for instance, to the discontinuous columns 



of damaged material, e.g., YBa2Cu3C>7 irradiated with 0.58 GeV Sn ions [13[], to carbon 



nanotubes embedded in Bi 2 Sr2CaCu 2 :c WM or to MgO nanorods grown in BSCCO super- 



conductors [15[|. The pinning force with a range of the order of £ is taken to be infinitely large 



- a trapped vortex segment cannot escape from the pinning center - and, therefore, only the 
interaction between the vortices is relevant. This assumption is physically reasonable since 
the maximal force the vortices exert on the pinned vortex segment is much smaller than the 
maximal possible force resulting from the depairing current. The driving force acting on 
the vortices is in the plane perpendicular to the VL and the magnetic field is restricted to 



B < 0.2H c2} i.e., the London approximation can be used, and to A > a = y^ /B, where A 
is the London penetration depth. These restrictions on the magnetic field are not severe as 
most of the experimentally accessible regime is covered. Within this regime, the scaling rules 
Tfifl are applicable allowing to generalize the results obtained for isotropic superconductors 



considered in the following. For convenience and technical reasons the moving VL is chosen 
as the frame of reference and hence the pinning center appears to move through the VL. 
In order to describe the lattice deformation produced by a defect we allow a finite number 
of soft vortices in the neighborhood of the pinned vortex line to accommodate according to 
their interaction. All other vortices are held fixed and consequently denoted as hard vortices. 
Accordingly, we split the (isotropic) London energy functional for the soft vortices 

P = i E/fri-dr A _ , (1) 



into three different parts: The self interaction accounting for the line energy of the soft 
vortices, the mutual interaction of the soft vortices, and the coupling of the soft vortices 
to the surrounding hard vortices. The instability to fluctuations of the energy functional 
([!]) is cured by taking a core term into account and by using a cutoff ~ £. The algorithm 
to relax the vortices is based on a conjugate gradient method; for details on the numerics, 



see Ref. |[L7| . We do not investigate the mechanism of vortex trapping and, therefore, the 
pinning center considered is initially positioned at a lattice site. The pinning center fixes 
the vortex along its length l p . Subsequently, the pinning center together with the attached 
vortex segment is adiabatically dragged through the VL while the deformation energy of the 
VL is monitored. The parameters in this procedure are the dragging angle d with respect 
to a basis lattice vector, the distance d from the equilibrium lattice site, and the length l p of 
the pinning center. The energy per length of the lattice deformations is measured in units 
of £o> where Sq = ($o/47rA) 2 , and the lattice spacing oq is used as the length unit. 

The point-like pinning centers have an extension of 3£ along the field direction and the 
magnetic field B pa 0.0025if C 2 is fixed by setting ao/£ = 50. Note that the pinning length 
l p = 3£ is larger than the cutoff £ of the energy functional but still small enough to correspond 
to point-like pins. Quantitatively similar results are obtained when l p is doubled. Estimates 
on elastic deformations of the VL suggest that the set of soft vortices should include the 
pinned vortex as well as the nearest neighboring vortices. Small displacements of the pinned 
vortex segment entail pronounced deformations only close to the pin where sharps kinks 
occur. With increasing displacement d the parts of the pinned vortex above and below the 
pin tend to align anti-parallel and, therefore, attract each other, see Fig. |I|. At the critical 
distance, dd p = 0.22ao, this attraction dominates and makes the configuration collapse, i.e., 
the two anti-parallel vortex segments annihilate, leaving behind a free unpinned vortex and a 
vortex loop attached to the pin. Subsequently, the free vortex line relaxes to its equilibrium 
lattice position and the vortex loop shrinks and disappears. With this one-loop depinning 
process the pin is freed and hence the VL becomes unpinned. Note that the change in 
topology is mediated by vortex cutting and recombination. The deformation energy versus 
the displacement is shown in Fig. [I| The sharp edge of the deformation energy at the 
recombination of the vortices, see Fig. [j], strongly influences the exponent in the scaling 
law of the depinning barrier and, in particular, seems to be the signature of depinning from 
strong pins. As the depinning distance dd p = 0.22ao is rather small, the pinned vortex line 
does neither noticeably affect the other vortices during the dragging process nor experience 



the non-circularity of the hexagonal VL potential. Consequently, the whole dragging process 
is essentially circular symmetric rendering the energy of the lattice deformation as well as 
the depinning mechanism independent of the angle. As a result, the presence of point-like 
pinning centers does not ensue a preferred orientation for a moving VL and, in addition, no 
vortex entanglement is introduced. 
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FIG. 1. The deformation energy as a function of the displacement of the pinned vortex segment 
is independent of the angle $ (the curves = 0°, 15°, and 30° coincide). The magnetic field is 
fixed at B ~ 0.0025i? C 2 (i.e., ao/£ = 50) and the pinning length is 3^. The depinning occurs at 
ddp = 0.22ao and the maximal pinning force amounts to F p ; n = 1.58£q. 

The critical current density j c derives from the critical state and, therefore, is related to 
the maximal pinning force per defect, F pin , by |18 



cnF, 



pin 



B 



(2) 



Here, n is the volume concentration of effective pinning centers which depends on the total 
defect concentration and the effective trapping area S according to n — UoS/Sq, where 
So = (V3/2)a 2 is the of the unit cell. On the basis of the depinning distance dd p an 
upper limit for S is obtained by 



s 




0) 



(3) 



So, otherwise. 

Inserting the numerical values F pin = 1.58e and d dp = 0.22a the upper limit for the critical 
current density becomes 



where j = (4/3^3) (ce )/(£$o) is the depairing current density. Note that j c is inversely 
proportional to B. For the creep activation energy, U(j), we obtain 



with U c = 0.18 ± 0. 02^0^0 an d a ~ 1.97. The exponent a ~ 2, resulting from the sharp 
drop in energy at depinning, differs from the value 5/4 obtained for a string in a (smooth) 
washboard pinning potential |i"9fl . 

Elastic considerations suggest that extended pinning centers deform not only the nearest 
neighboring vortices but also some of the next nearest neighbors and, therefore, all these 
vortices must be soft. In order to study the effect of extended pins, the magnetic field is 
fixed at B w 0.06i/ C 2 corresponding to ao/£ = 10 and the pinning centers are restricted to 
lengths l p < 7a Q . Firstly, we analyze the topological aspects of the lattice deformations and, 
secondly, we determine the critical current density and the activation energy for depinning. 

For extended pinning centers the situation is more complicated than for point pins and 
the lattice deformations produced depend crucially on the angle the pinning length l p , 
and the distance d. In order to illustrate the effects of extended pins we consider a pin with 
a length l v = 3ao dragged along an angle $ = 5°. At a distance d = 0.98a configurations 
of four twisted vortices, so-called twisted quadruplets (TQ), are formed above and below the 
pinning center, i.e., a loop and anti-loop excitation of four vortices, see Fig. 0. 
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FIG. 2. An extended pinning center with l p = 3ao together with the pinned vortex segment is 
dragged through the VL at an angle # = 5°. Two TQ's appear at = 0.98ao and disappear again 
after the vortex has depinned via a pair annihilation state. 

At the transition from the disentangled to the entangled state the deformation energy drops 
considerably ~ (0.5 — l)e a , see Fig. [5], and becomes slightly larger than the energy of 
two TQs because the pin is not at a lattice site. The high energy of these loop excitations, 
E T q — E = 2Ale a (E is the ground state energy of the regular VL), is due to the 
hexagonal VL which confines the TQ to a length 1?q = 2.94ao- The TQs above and below the 
pinning center have opposite orientation and, therefore, attract each other. The interaction 
energy of the TQs can be estimated by considering two oppositely oriented vortex rings in 
parallel planes separated by a distance L. Choosing vortex rings on top of each other with 
the same radius R, we obtain 




(6) 



Using the numerically obtained values L m 6.4a and R ~ a we find 



E L (TQ) « -2n 2 e a A /L 3 « -0.08e ao- (7) 

When the pinned vortex segment is at its closest distance from its nearest neighboring vortex, 
the configuration is similar to the initial state since the TQs "carry away" the displacement 
of one lattice constant. Therefore, it looks as if the nearest neighbor rather than the original 
vortex were pinned. 
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FIG. 3. The deformation energy exhibits a palpable drop when the topology changes, i.e., for 
a displacement d ~ ao- For $ = 0° a one-loop depinning occurs whereas for small angles, e.g., 
$ = 5°, TQs are generated. Intermediate angles, e.g., $ = 8° and •& = 15° lead to the creation of 
TTs. For larger angles i? = 18° and $ = 28° a pair annihilation state is formed. Note that the 
largest slope of the energy curve, obtained for # = 0°, determines the maximal pinning force and 
hence the critical current density as well as the orientation of motion. The configurations (a), (6), 
and (c) are displayed in Fig. ||. 

If the pinned vortex segment is dragged beyond the nearest neighboring vortex towards 
the next nearest neighbor the distance between the TQs is increased. Despite the larger 



S 



separation of the TQs, there is not enough space for the creation of a second pair of TQs (this 
behavior has been obtained for l p < 7a ). Instead, the pinned vortex and the next nearest 
neighboring vortex form a configuration of two twisted vortices. Since configurations of two 



twisted vortices are unstable in the VL |T7| the two vortices collapse and a pair annihilation 
state is formed, leading to an exchange of the pinned vortex segment, see Fig. After 
depinning the two oppositely oriented TQs move together and annihilate, leaving behind 
a "healed" regular VL. The depinning terminates the generation of vortex entanglement 
associated with the formerly pinned vortex. However, configurations of entangled vortices 
can be stabilized by other pins, in particular, weak point disorder. 

The formation of two TQs is a generic case for extended pinning centers. With increasing 
angle, the force the pinned vortex exerts onto its closest neighboring vortex is reduced and, 
therefore, only configurations of three twisted vortices, so-called twisted triplets (TT), are 
generated. The TTs behave similarly to the TQs and, in particular, disappear in the same 
way. For angles close to tt/Q a pair annihilation state is formed (note that the twisted pair is 
unstable in the VL) which in turn leads to an exchange of the pinned vortex segment with 
the nearest neighboring vortex. The lattice deformations as a function of the angle $ and 
the pinning length l p are summarized in the "phase diagram" in Fig. 
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FIG. 4. Extended pinning centers can trigger the formation of configurations with three or four 
twisted vortices. Below the solid line dep inning is mediated by a vortex loop creation at the pin, 
whereas above this line depinning always involves pair annihilation either with the nearest or the 
next nearest neighboring vortex. 

The strong dependence of the deformation energy on the angle see Fig. ^ makes the 
VL move along a lattice vector. This behavior is a consequence of the maximal pinning force 
which is obtained for $ = 0°, see Fig. |3|. Interestingly, weak pinning seems to introduce the 

as 



same orientation of motion as has been demonstrated by experiments on Al films [20 



well as by theoretical studies [21,22]. The critical current density for extended pins derives 
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along the same lines as for point pins, but with S = So, as the change in topology occurs at 
d ~ a . Accordingly, the critical current density is given by 

3^/3 F p - m 2 . 

Jc = — A "oOoOo- ( 8 ) 

F pin depends on the field and, therefore, the critical current density is not inversely propor- 
tional to B (as is the case for point pins) but is given by 

iV 7 



j c «6.l(^J n ^jo, (9) 

with 7 = 0.7 ± 0.1, where B = $o/«o * s normalized with respect to B ip = $oAp- The 
activation energy to change the topology, however, exhibits the same scaling behavior as for 
point-pins 

uu) = um(^-J , (io) 

with a « 2, U c (l p ) = f3e a (l p /a )~< , f3 = 1.63 ± 0.05 and 7 = 0.78 ± 0.05. The universality 
of the exponent a originates from the pronounced edge of the deformation energy when the 
topology changes. 

In conclusion, we have studied the plastic deformations of the vortex lines due to strong 
pinning. We have identified the relevant depinning processes involving either loop creation 
or pair-annihilation, depending on the pin size. We have found that large pins do generate 
entangled configurations (TTs and TQs) which will be stabilized by a finite density of pins. 
Finally, we have determined the critical current density j c as well as the activation barrier 
U(j) for creep, the latter showing a universal exponent a ~ 2 originating from the change 
in topology at depinning. 
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